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ABSTRACT 

We examine standard methods of measuring nebular chemical abundances, including estimates based 
on direct T e measurements, and also emission-line diagnostics. We use observations of the LMC H n 
regions DEM L199, DEM L243, DEM L301, and DEM L323, whose ionizing stars have classifications 
ranging from 07 to WN3. Following common practice, we assume a two-zone T e structure given by 
T(0 ++ ) and T(0 + ) to compute ionic abundances. We compare with photoionization models tailored to 
the observed properties of the individual objects, and emphasize the importance of correctly relating T e 
in the two zones, which can otherwise cause errors of ~0.2 dex in abundance estimates. The data show 
no spatial variations or local metallicity enhancements to within 0.1 - 0.15 dex in any of the objects, 
notably including DEM L199, which hosts three WR stars. 

Our data agree well with both the modeled i?23 and 523 abundance diagnostics for O and S. We present 
the first theoretical tracks for 523, which are in excellent agreement with a larger available dataset. 
However, contrary to earlier suggestions, 5*23 is much more sensitive to the ionization parameter (U) 
than is i?23. This occurs because 523 does not sample S IV, which is often a significant population. We 
therefore introduce 5234 =([S II] + [S III] + [S IV])/H/3, and demonstrate that it is virtually independent 
of U. Predicted and observed spatial variations in 5234 are thus dramatically decreased in contrast to 
523. The intensity of [S IV]A10.5^ can be easily estimated from the simple correspondence between 
[S IV] /[S III] and [O III]/[0 II]. Using this method to estimate 5234 for data in the literature yields 
excellent agreement with our model tracks, hence we give a theoretical calibration for 5234. Our models 
show that the double-valued structure of 523 and 5234 remains an important problem as for i?23, and 
presently we consider calibrations of these S diagnostics reliable only at Z <; 0.5 Zq. However, the slightly 
larger dynamic range and excellent compatibility with theoretical predictions suggest the S parameters 
to be more effective abundance diagnostics than i?23. 

Subject headings: galaxies: abundances — galaxies: ISM — Hn regions — ISM: abundances — 
Magellanic Clouds — supernova remnants 



1. INTRODUCTION 

The emission-line signatures of H II regions are a power- 
ful and widely-used indicator of galactic abundances. This 
is especially true in distant galaxies where most stellar 
abundance probes cannot be employed. Furthermore, the 
spectral properties of H n regions also give important diag- 
nostics of the ionizing stellar population, such as effective 
temperature and numbers of stars. Given the wide use of 
such nebular diagnostics, it is vital to test and calibrate 
them using H n regions with known characteristics and 
ionizing stellar populations. 

Along with a companion paper (Oey etal. 2000; here- 
after Paper I), we report here on a detailed study for this 
purpose, of four H n regions in the Large Magellanic Cloud 
(LMC). The OB associations in all four of these Hn re- 
gions have been examined in detail and classified, thereby 
strongly constraining the ionizing energy distributions. In 
addition, the LMC's proximity and orientation with re- 
spect to the Galaxy also permit a detailed understanding 
of the nebular morphology. Thirdly, the abundances can 
be accurately determined. We therefore have high-quality 
information on the three principal parameters that deter- 
mine the nebular emission: stellar effective temperature 



(T*), the ionization parameter (U) that relates ionizing 
flux to gas density, and metallicity (Z). 

Table m gives a brief summary of our objects, which can 
be examined in greater detail in Paper I. The first two 
columns in Table [l] give the H n region identification in 
the Davies, Elliott, & Meaburn (1976) and Henize (1956) 
Ha catalogs, respectively; the third column identifies the 
parent OB association from Lucke & Hodge (1970). Col- 
umn four shows the spectral type of the dominant ionizing 
stars, as classified by the references in column 8. Column 
five lists the nebular Ha luminosity from Oey & Kenni- 
cutt (1997), and column 6 shows the adopted value from 
Paper I, of the inner radius of the gas distribution, as a 
fraction of the Stromgren radius Rs- This gives some in- 
dication of the nebular morphology. Finally, column 7 in- 
dicates the presence of significant shock excitation: DEM 
L243 includes an embedded or superimposed supernova 
remnant (SNR), and DEM L301 also shows evidence of 
significant shock activity (Paper I). While shocks may be 
present in the other objects as well, their contribution to 
the nebular emission is unimportant. 

The detailed presentation of the objects is given in Pa- 
per I, including narrow-band images in Ha, [O III], and 
[S II] , and spectroscopic observations over the wavelength 
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range 3500 - 9200 A. For each object, we observed two to 
three stationary, spatially-resolved slit positions. For all of 
the objects except DEM L301, we also obtained spatially 
integrated observations by scanning the long slit across 
the nebulae. The scanned data should therefore resemble 
typical observations of such H n regions at distances of 10 
- 20 Mpc. As seen in Table [l], the sample spans a range 
in dominant stellar spectral types, from 07 to early Wolf- 
Rayet (WR). There is also variety in the morphology of 
the objects, ranging from near perfect Stromgren sphere 
(DEM L323), to extreme shell structure (DEM L301), to 
highly complex (DEM L199). 

Paper I provides a detailed analysis of the spatially re- 
solved emission-line ratios with respect to the sequence in 
represented by these objects, and a comparison with 
photoionization models using the current generation of 
stellar atmosphere models for both O stars and early WR 
stars. In general we found a gratifyingly high level of agree- 
ment, largely supporting the CoStar energy distributions 
of Schaerer & de Koter (1997) and early WR atmospheres 
of Schmutz, Leitherer, & Gruenwald (1992) and Hamann 
& Koesterke (1998). In conjunction with data from Kenni- 
cutt etal. (2000), we provide a first, empirical calibration 
of nebular diagnostics for the dominant T+. In addition to 
the well-known rf radiation softness parameter of Vflchcz 
& Pagel (1988), we introduce [Ne IIl]A3869/H/3 as a simi- 
lar diagnostic, which is more robust to nebular conditions 
and is sensitive for a higher range of T* . 

We also presented in Paper I the spatially-resolved be- 
havior of the O abundance parameter (Pagel etal. 1979), 



i?23 



_ [O II]A3727 + [O III]AA4959, 5007 



(1) 



and S abundance parameter (Vflchez & Esteban 1996; 
Christensen etal. 1997; Diaz etal. 1999), 



523 



_ [SII]A6724+ [S III] AA9069, 9532 
= H/3 



(2) 



where we designate the [S II] lines A6716 + A6732 as A6724, 
analogous to [O II]A3727. We confirmed that both obser- 
vationally and theoretically, R23 remains spatially uniform 
across the nebulae. In contrast, for uniform abundances, 
models of 523 vary across the nebulae, showing lower val- 
ues in central regions. The observations clearly reflect this 
pattern, which is caused by the missing ionization stage of 
S IV. We address this issue in detail below in this paper. 
As is well-known, i?23 and 523 are also sensitive to T*, 
and these effects are also shown in Paper I for our objects. 

In this paper, we present the abundance determinations 
for our sample, for both the spatially-resolved and scanned 
longslit observations. We examine conventional assump- 
tions for the nebular electron temperature (T e ) structure 
and resulting ionic and total abundance determinations. 
We then explore the metallicity diagnostics i?23 and 523 
in more detail, and introduce another diagnostic, 5234. 
As before, we use the photoionization code MAPPINGS II 
(Sutherland & Dopita 1993) in conjunction with CoStar 
stellar atmosphere models (Schaerer & de Koter 1997) for 
O stars. 



2. DIRECT ABUNDANCE DETERMINATIONS 

We derived abundances from the measured line fluxes 
in Paper I using standard techniques. For calculation of 
line emissivities for elements heavier than helium, we used 
the five-level atom code Five_L as implemented in Stsdas 
version 1.3.5 (Shaw & Dufour 1995). We initially obtained 
an estimate of electron density from the [S II] A6716/A6731 
line intensity ratio, assuming a default T e = 10 4 K. The 
densities are almost all ^ 100 cm~ 3 . 

2.1. Temperature structure 

The assumed T e structure, and thereby, ionization struc- 
ture of the H II region, can generate substantial uncer- 
tainties in even "direct" abundance determinations from 
typical optical emission-line spectra. The T e structure is 
in turn a function of density structure and stellar atmo- 
sphere models. Peimbert's (1967) t 2 formulation for tem- 
perature fluctuations is probably the best-known example 
of addressing this problem, which can be used for both 
small-scale and large-scale variations in T e . 

We consider here the large-scale temperature structure. 
As discussed below, our photoionization models tend to 
overestimate T e , suggesting that small-scale fluctuations 
are not significant. Standard abundance determinations 
adopt either a single characteristic T e or a two-zone model 
for T e , and these assumptions can also cause substantial 
errors in abundance determinations. Garnett (1992) pro- 
vides one of the more thorough investigations of this issue, 
which is a well-known problem for elements like S that 
do not conform well to a two-zone model. However, if 
the adopted zone temperatures do not adequately charac- 
terize the respective regions, then the inferred abundance 
can be significantly in error even for elements like O that 
are well-described by a two-zone model. This can be a 
problem especially in cooler nebulae with T e <; 10, 000 
K, that have strong temperature gradients (Garnett 1992; 
Stasihska 1980). 

We assumed here a two-zone temperature structure for 
our nebulae, such that a common T e was adopted for ex- 
citation of O III, Ne III, S III, and Ar III; while a separate 
T e was adopted for the excitation of N II, O II, and S II. 
We used the same electron density for both zones. 

For the high ionization zone, we adopted T(0 ++ ), the 
volume-averaged T e for the ++ population. This is ob- 
tained essentially directly from the observed T[0 III], the 
temperature inferred from the [O III] A4363/(A4959 + 
A5007) line intensity ratio. A temperature can also be ob- 
tained in principle from the [S III] A6312/(A9069 + A9532) 
ratio. However, the auroral line flux in this case was of- 
ten described by a large fractional uncertainty, which in- 
troduced correspondingly large random errors in the re- 
sulting abundances; there is also some possibility that 
the near-IR lines may be affected by telluric absorption 
(e.g., Stevenson 1994). We consequently chose to adopt 
T(0 ++ )~ T[0 III] for the high-ionization zone, includ- 
ing the S ++ region. Garnett (1992) has suggested on the 
basis of nebular models that the relation between ion- 
weighted temperatures for S ++ and ++ is linear, but 
with a slope differing from unity. For objects in the current 
study, T[0 III] values are generally close to 10,000 K, at 
which point Garnett 's work suggests that the ion- weighted 
T(S++) and T(0++) should be similar; thermal effects of 
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dust in the nebula may also be expected to reduce the con- 
trast between the two temperatures (Shields & Kennicutt 
1995), lending support to use of a common value. 

For the low-ionization zone, T e Ccin b6 obtained directly 
from the [S II](A4069 + A4076)/(A6716 + A6731) ratio, 
but we found that this option frequently suffered from 
a low signal-to-noise ratio. It is often standard prac- 
tice to derive the lower-ionization T(0 + ) from the higher- 
ionization T(0 ++ ) using an analytic relation (e.g., Camp- 
bell, Terlevich, & Melnick 1986; Pagel etal. 1992), whose 
accuracy clearly affects that of the abundance determina- 
tion. We investigated several alternatives before adopting 
a prescription for T(0 + ). In deciding what relation to 
adopt, our principal criterion was that the abundances in- 
put to Mappings photoionization models for the individ- 
ual objects (Paper I) should be consistent with those in- 
ferred from the output emission-line spectra. For example, 
when using the expression employed by Pagel et al. (1992, 
their equation 6), our resultant Hn region model spec- 
tra implied an O abundance that was ~ 0.2 dex lower 
than the assumed input, a substantial departure from self- 
consistency! 

The relationship between T(0 ++ ) and T(0 + ) is model- 
dependent, and varies with ionization parameter and 
metallicity. Figure shows T(0+) vs. T(0++) for Map- 
pings photoionization models at Z = 0.3 Z©, with an in- 
ner radius of 0.4i?s- The different symbols correspond 
to different CoStar stellar atmospheres as indicated, with 
B2, C2, and E2 corresponding to spectral types 08-09, 
06-07, and 03-04, respectively. We show models with 
\ogU = —2, —3, and —4. The dashed and dot-dashed 
lines show the relation between T(0 + ) and T(0 ++ ) from 
Campbell etal. (1986; see also Garnett 1992), and Pagel 
etal. (1992), respectively, while the dotted line delineates 
T(0+)=T(0++). 

The Pagel et al. relation deviates significantly from the 
models at these temperatures, while the Campbell et al. re- 
lation shows reasonable agreement in slope for T(0 ++ ) 
> 10,000 K, falling close to the theoretical predictions for 
logt/ = —3. At lower temperatures, the models are more 
consistent with an isothermal nebula. We consequently 
adopted the formulation, 



T(0+) = / °' 7 T(0+ ~ 



3000 K 



T(0 H 



T(0++) > 10,000 K 
T(0++) < 10,000 K 

(3) 

where T(0++)= T[0 III] as described above. The re- 
lation for T(0++) > 10,000 K is that given by Camp- 
bell et al. Equation |^ yields consistent input and output 
abundances to < 0.05 dex, except for Ne, whose modeled 
line ratios are persistently discrepant with the observa- 
tions (Paper I; see below). We caution that Figure [j] and 
equation are optimized in the parameter space relevant 
to our LMC targets. Different formulations may be more 
appropriate at other metallicities and ionization parame- 
ters; it is beyond the scope of this work to fully examine 
this issue. However, it is clear that some care is necessary 
in choosing a relation between T(0 ++ ) and T(0 + ). 

While we ensured consistency between the input and 
output abundances to the photoionization models, one 
worrisome problem that remains unresolved is the over- 
prediction of T e in comparison to the observations, as pre- 
sented in Paper I. To briefly summarize, the predicted 



T[0 III] is greater than that inferred from the observed 
[O III]A4363/A5007 by 850 K in DEM L323 and 1500 K 
in DEM L199. The problem may also exist in the other 
two objects, DEM L243 and DEM L301, but the fainter 
emission from A4363 in these lower-excitation nebulae pre- 
vented any useful constraints. We also found that the dis- 
crepancy persisted when using Hummer & Mihalas (1970) 
stellar atmosphere models and the Cloudy (Ferland 1998) 
photoionization code. Since T(0 ++ ) and T[0 III] agree 
to < 1% in the models, the effect is not caused by non- 
collisional excitation of A4363 in the models. As discussed 
in Paper I, the temperature discrepancy occurs in lines 
of sight across the entire nebulae, and is in the opposite 
sense of that expected from small-scale temperature fluc- 
tuations (Peimbert 1967). Our sample therefore shows no 
evidence that such T e fluctuations systematically bias our 
abundance determinations. Mathis (1995) discusses evi- 
dence for and against the general existence of significant 
T e fluctuations. 

2.2. Ionic abundances 

Determining ionic abundances relative to H requires 
calculation of the H/3 emissivity, which is also tempera- 
ture dependent, but originates in both the high- and low- 
ionization zones. We adopted an intermediate tempera- 
ture for calculation of the H/3 emissivity, representing an 
average of T(0 ++ ) and T(0 + ), weighted by the relative 
abundances of ++ and + . The same intermediate tem- 
perature was used for computing He emissivities, using 
expressions taken from Benjamin et al. (1999). This con- 
struction was used to obtain ionic abundances from the 
line measurements in Paper I, which were corrected for 
reddening and, where necessary, underlying Balmer ab- 
sorption. The abundance results for different nebular lines 
of a single ion were combined, weighted by the variance of 
the line fluxes. 

Total abundances were obtained from the ionic abun- 
dances following the standard practice of using ionization 
correction factors (ICFs) to allow for unobserved ioniza- 
tion stages. We employ the following relations: 



O _ 0+ + Q++ 
H ~~ H+ 



(4) 



where the lack of detectable He IIA4686 emission indicates 
no significant +3 population; 



(Peimbert 



N _ N+ O 
H 0+ H 
Costero 1969; Garnett 1990); and 

~S+ + S+ 



where 



ICF(S+ 3 ) 



1 



'lCF(S +d ) 



1 



0+ 

o 



1/c 



(5) 



(6) 



(7) 



with a = 2.5, corrects for unobserved S +3 ions (Garnett 
1989; Stasihska 1978). For the noble gases, we adopt 
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Fig. 1. — T(0+) vs. T(0 ++ ) from photoionization models. The symbols correspond to models using the indicated CoStar stellar atmo- 
spheres (see text). The locus of results for log U = —2, —3, and —4 fall in the sense indicated for the C2 runs. We overplot dashed and 
dot-dashed lines that show the relations of Campbell et al. (1986) and Pagel et al. (1992), respectively; the dotted line shows T(0 + )=T(0 ++ ). 



(Pcimbert & Costero 1969; Simpson etal. 1995), and 

Ar _ Ar++ S 
IT ~ S+ + S++ ' H 



(9) 



(Garnett etal. 1997). Most of our objects are ionized by 
early-type stars (Table lj), thereby fully ionizing He, but 
not exhibiting detectable He II. Thus our default relation 
for He is simply, 

He He+ 

TT-1F ' (10) 

For DEM L243, which is ionized by 07 stars, the He I 
A5876/H/3 line ratios suggest that He is not fully ionized 
throughout the nebula (Paper I; see below). For this ob- 
ject, we experimented with the expression for He/H given 
by Peimbert & Torres-Peimbert (1977; their equation 15), 
that includes an ICF for He ; however, this prescription 
yields He abundances that are at least 0.1 dex higher than 
the rest of the objects in the sample. We therefore give re- 
sults for DEM L243 using equation [lO] and note that these 
should be treated as lower limits to the true abundance. 
It is worth remarking that our investigation supports the 
results of Baldwin etal. (1991), who find that the Orion 
He abundance was slightly overestimated by Peimbert & 
Torres-Peimbert as a result of their ICF for He . 

The theoretical support for use of equations || - ^ derives 
mostly from models that consider the integrated prop- 
erties of Hn regions (e.g., Garnett 1992; Mathis 1985; 
Stasinska 1978). Thus, observations that sample only a 
small "pencil-beam" through a nebula may not yield the 
correct abundances if such an analysis scheme is used; Gru- 
enwald & Viegas (1992) have discussed this problem in de- 



tail. However, previous observational studies have found 
generally good agreement in elemental abundances de- 
rived at variable positions across individual nebulae, using 
the integrated-spectrum methods (e.g., Diaz etal. 1987; 
Masegosa, Moles, & del Olmo 1991; Gonzalez-Delgado 
etal. 1994). Such agreement was found even when sub- 
stantial excitation gradients were seen within the H n re- 
gion. In the present study, we can directly compare abun- 
dances obtained from the integrated spectrum of a neb- 
ula with those calculated by the same means from small- 
aperture measurements at a variety of radii. We show 
below in §2.3.1 that we find a high degree of consistency 
between the different aperture measurements for the same 
source, including the integrated spectrum. These results 
support the validity of using integrated-spectrum prescrip- 
tions when only a part of the H n region is observed. Ap- 
parently the characteristic physical parameters determined 
for the two ionization zones are adequate to determine the 
ionic abundances to high accuracy in our objects. 

2.3. Abundance results 

We used the prescriptions above to compute abundances 
for both the spatially resolved and scanned, spatially in- 
tegrated observations for each object. Table || presents 
results for the individual, spatially resolved observations, 
along with error estimates. These were obtained from 
the uncertainties in measured fluxes, using a Monte Carlo 
method similar to that described by Kobulnicky & Skill- 
man (1996). Input flux values were modified randomly 
by addition of Gaussian noise with amplitude correspond- 
ing to 1-er measurement uncertainties, for a total of 5000 
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iterations per spectrum. A few points with low signal- 
to-noise ratios yielded unphysical solutions in the T e and 
abundance determinations, and these were discarded from 
the sample. The standard deviation of the resulting abun- 
dance distribution was adopted as the final uncertainty. 

Table || presents the mean values of the spatially- 
resolved measurements, weighted inversely by the vari- 
ances. We also list the corresponding formal uncertainty 
in the mean. Note that the weighting scheme and quoted 
errors assume that actual variations are negligible in the 
quantities listed in Table ||, and that the scatter results 
from a normal distribution of measurement errors. As 
described below, we do not find significant evidence for 
abundance fluctuations; however, it is possible that the 
measurement errors may not be a strictly normal distri- 
bution. The listed errors in Table || are therefore likely 
to somewhat underestimate the true uncertainties, nor do 
these include systematic errors. 

Regarding the abundances of Ne, it is important to 
note that there is probably a substantial systematic er- 
ror in the values obtained in Table ^, and also later in 
Table ||. As discussed in Paper I, the observed line intensi- 
ties for Ne are systematically discrepant with the tailored 
models. Although earlier generation stellar atmospheres 
caused [Ne III] intensities to be underpredicted, we now 
find a modest overprediction. The discrepancy can be re- 
solved by reducing log(Ne/H) by 0.2 - 0.3 dex. However, 
since we do not understand the cause of the discrepancy, 
we have chosen to list the Ne abundance derived from the 
standard relations. It appears likely that errors in the 
stellar atmospheres are responsible for much of the prob- 
lem (Paper I). Peimbert (1993) also emphasizes the uncer- 
tainty in equation ||. The high ionization potential (40.96 
eV) required for Ne III gives it outstanding potential for 
probing hot stellar ionizing sources, so it is highly desire- 
able to resolve the uncertainties regarding its emission. 

2.3.1. Spatial uniformity 

In Figure ^p, — e, we show the spatial distribution of ele- 
mental abundances across the sightlines for our stationary 
observations of DEM L199. Figure ||/ shows the corre- 
sponding distribution in T[0 III] measurements. For DEM 
L243, DEM L301, and DEM L323, we show in Figures |~| 
the results for He, N, and O, as representative of the other 
elements. These are shown respectively in panels a, b, 
and c; and panels d show corresponding measurements of 
T[0 III]. The different symbols correspond to individual 
slit positions as designated in Paper I (Figures 6, 8, 12, 
and 13), to aid cross-referencing. Following our conven- 
tion in that work, we simply show the slit positions su- 
perposed, therefore these figures will not necessarily show 
clean radial profiles across the nebulae, although they do 
approximate this reasonably well. This issue may be in- 
spected in Paper I, along with actual slit positions. The 
light, horizontal lines indicate the spatial extent of each 
extracted aperture, while the vertical error bars indicate 
uncertainties in the derived measurements. 

In Figure |L DEM L243 shows much larger scatter and 
uncertainty among the individual apertures, which can 
be attributed to the large measurement uncertainties for 
T[0 III] (Figure [§f). There is also an SNR either em- 
bedded or superimposed in the line of sight to this ob- 
ject, but we have not included the affected apertures in 



our abundance estimates or in Figure g[ Assuming an 
essentially constant metallicity in DEM L243, we recom- 
puted the abundances by fixing T(0 ++ ) to be the mean 
T[0 III], using the same weighting and omitting the same 
deviant values as before, in computing the mean abun- 
dances. Figure ^ and Table || show the results using the 
new T(O ++ )=9700 K. The dramatic reduction in scatter 
for the heavy elements suggests that the original scatter 
in Figure may indeed be caused by poor measurements 
in T[0 III] . However, we emphasize a point by Mathis, 
Chu, & Peterson (1985), that true abundance fluctuations 
will induce corresponding fluctuations in T e , owing to the 
more efficient cooling accompanying higher O/H. The re- 
duced scatter with a fixed T e is therefore not a strong 
demonstration of truly constant abundances. To further 
test consistency with constant abundance distributions, we 
computed the reduced-^ 2 statistics for the original distri- 
bution, yielding 1.4 and 1.8 for log(N/H) and log(0/H), 
respectively. The probabilities of obtaining these values 
for 12 degrees of freedom are 16% and 4%, respectively. 
These values therefore hint that the scatter may in part 
be caused by real variations, although the significance is 
low. 

In the event that the abundances for DEM L243 are 
essentially constant, the mean values obtained with fixed 
T[0 III] should give a better estimate than the original 
values (Table |). However, we caution that there is sys- 
tematic uncertainty introduced by the adopted value of 
T(0 + +); comparison to the original mean abundances in 
Table || suggests the uncertainty is roughly 0.1 dex in the 
derived metallicities. Since the recombination lines used to 
determine log(He/H) are less temperature-sensitive than 
the collisional metal lines, log(He/H) (Figure ||a) still 
shows larger scatter than exhibited in the other objects. 
However, note that the magnitude of the variation is only 
~ 0.3 dex. As discussed in Paper I, the He I A5876/H/3 ra- 
tios suggest that He is not uniformly fully ionized in DEM 
L243, therefore causing log(He/H) to be underestimated in 
many apertures. Our value for the He abundance is there- 
fore a lower limit in this object. The upper envelope to 
the distribution in Figure ||a is around log(He/H) ~ —1.1, 
which may be more indicative of the true He abundance. 
This value is more consistent with those for the other ob- 
jects (Table |). 

Thus, in Figures |^ the abundances appear to show no 
spatial variations within the measurement uncertainties, 
with standard deviations typically around 0.10 - 0.15 dex. 
The uniformity of the abundance derivations is reassuring, 
and suggests that our adopted ionic relations (§2.2) and 
description of the nebular temperature structure (equa- 
tion ||) yield abundance estimates with high accuracy. The 
apparent success of these methods is consistent with the 
finding by Mathis (1985) that ICFs for model nebulae ap- 
pear fairly robust between volume averaged and radially 
averaged regions. The uniformity of results for apertures 
tracing different parts of the same nebula further suggests 
that spurious "pencil-beam" effects resulting from projec- 
tion of radial T e gradients (Gruenwald & Viegas 1992) 
are negligible for our objects. This is consistent with the 
fact that sources in our sample have T[0 III] restricted 
to the approximate range 9500 - 12,000 K, for which the 
nebulae are expected to be relatively isothermal (equa- 
tion ||; Figure |l|). We emphasize that this propitious cir- 
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Fig. 2. — Spatial variation in abundances and T[0 III] along the slit positions for DEM L199. Light, horizontal lines show aperture sizes, 
while vertical error bars indicate abundance uncertainties. An LMC distance of 50 kpc corresponds to 0.24 pc arcsec -1 . Apertures near WR 
stars are highlighted with solid dots. 
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Fig. 6.— Data for DEM L243 as in Figure Bl but with T(0++) fixed at 9700 K for all included apertures. 
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cumstance may not apply to other Hn regions, that may 
exhibit stronger T e gradients; as an example, Walter, Du- 
four, & Hester (1992) reported a significant T e gradient 
in the Orion nebula, although their spatial scale of 0.5 pc 
is almost an order of magnitude smaller than is relevant 
for our objects in the present study. As in §2.1 above, we 
again emphasize the importance of choosing the correct 
parameterization for the T e structure. 

2.3.2. Search for self- enrichment 

It is especially interesting to examine DEM LI 99 for 
possible abundance variations that are introduced by the 
three WR stars. Two of these are binaries with a WN3 
component (Breysacher 1981) and one is a WC4 + 06 
V-III binary (Moffat etal. 1990). The exact location of 
these stars with respect to the nebular gas distribution 
and slit positions may be examined in Paper I. We em- 
phasize that DEM L199 is not a WR ejecta shell, but is 
a large, luminous H II complex whose dominant ionizing 
stars are early WR stars. Walsh & Roy (1989) and Kob- 
ulnicky et al. (1997) identified two regions in the starburst 
galaxy NGC 5253 that appear to show enhanced N, which 
is suggested to result from self-enrichment by WR stars. 
Kobulnicky etal. (1997) point out that accompanying He 
enrichment would be expected from the WR sources, and 
acknowledge that the lack of He excess in these objects is 
puzzling. 

DEM L199 has three WR stars within a large central 
cavity in a luminous H n region. It would seem likely 
that, if self-enrichment from WR stars can be seen, it 
should be detected in our observations. In Figure ||, we 
indicate six apertures that are closest to the stars Br 32 
(WC4) and Br 33 (WN3), with solid dots. These aper- 
tures are identified in Paper I as: D199. 205-15, 16, and 
17; and D199.205N120-17, 18, and 19. The characteristic 
physical distances from the WR stars are about 5-15 
pc. Figure || shows that these points do not show any sign 
of abundance enhancements. Results from the full set of 
aperture measurements for DEM LI 99 are consistent with 
uniform abundances at the 90% confidence level, as indi- 
cated by reduced-x 2 values of 0.59 and 0.60, for log(N/H) 
and log(0/H) respectively, with 17 degrees of freedom. For 
He, in particular, the variation is constrained to be <; 0.1 
dex. Although abundance enhancements are seen in WR 
ejecta nebulae, our data suggest that self-enrichment by 
WR stars can be an extremely subtle phenomenon. It is 
possible that the processed material is heated to coronal 
temperatures within the superbubble and is thus optically 
undetectable. Alternatively, the WR phase in these stars 
may not yet have lasted long enough to produce significant 
enrichment in the surrounding environment. 

We were also interested to see whether the SNR in DEM 
L243 showed evidence of abundance anomalies. Abun- 
dance estimates from SNRs are necessarily more uncer- 
tain, owing to the more complicated radiative transfer, but 
meaningful estimates have been made by, e.g., Russell & 
Dopita (1990). These are based on matching the emission- 
line spectra with a shock code. In Paper I, we were able to 
find excellent agreement between MAPPINGS shock models 
and the observed emission from the shock-affected aper- 
tures, using the abundances derived from the uncontam- 
inated regions of the nebula listed in Table |. The SNR 
therefore shows no evidence of abundance anomalies, al- 



though strong emission from photoionized gas in the same 
line of sight prevents strong constraints. 

Likewise, the superbubble DEM L301 shows strong ev- 
idence of a recent SNR impact (Chu & Mac Low 1990; 
Oey 1996b; Paper I), and thus could conceivably show 
enrichment by massive star winds and supernova ejecta. 
However, its elemental abundances are often lower than 
for the other objects in the sample (Table ||). We cau- 
tion that shock-excitation can significantly affect the abun- 
dance determinations for this object (Paper I). Peimbert, 
Sarmiento, & Fierro (1991) showed that contamination by 
shock activity can cause abundances to be underestimated, 
especially for higher-ionization species. Our measurements 
for DEM L301 are consistent with this behavior, although 
in §2.3.3 we find that the metallicities for DEM L243 do 
not appear strongly affected by SNR contamination. 

2.3.3. Spatially integrated abundances 
and LMC metallicities 

We also derived abundances for the scanned, spatially 
integrated apertures using the same methods. These are 
presented in Table [|, along with mean LMC Hn region 
abundances from compilations by Dufour (1984). A more 
recent compilation by Garnett (1999) shows essentially the 
same values. Within the uncertainties, the mean abun- 
dances for each object from Table || agree with the deter- 
minations from the spatially integrated observations (Ta- 
ble |), although the offsets appear to be systematic across 
all elements for each object. This again points to uncer- 
tainties in the T(0 ++ ) determination, which can result 
from simple measurement errors, or factors related to the 
spatial integration of the line emission. We regard the 
values obtained from the mean of the spatially resolved 
data (Table |) to be more reliable than the single observa- 
tions from the integrated data. However, we caution that 
the ionic relations in §2.2 are based on spatially integrated 
models and observations, and this could introduce system- 
atic variations between the spatially resolved and scanned 
data. But taking the derived values and uncertainties at 
face value, the mean abundances of the stationary aper- 
tures should be somewhat more reliable. 

We include in Table | the abundances for DEM L243 
derived from both the total, scanned observation and the 
scanned observation with the SNR-contaminated region 
subtracted. Interestingly, there is no significant difference 
between these, although the data including the SNR do 
show the expected decrease in computed abundances (Pe- 
imbert etal. 1991). Thus, while DEM L301 showed sus- 
piciously low metallicity measurements attributable to ef- 
fects of shock emission, DEM L243 is an example where 
the SNR is not a significant factor. In Paper I, we also 
found that the two shock-affected objects exhibit different 
behavior in their line emission with respect to the pho- 
toionized regions, thereby demonstrating how shocks con- 
tribute in different ways to the spectra of host H n regions, 
depending on shock velocity and environment. 

We find a tendency for our measurements to be ~ 0.2 
dex lower than the mean LMC metal abundances com- 
piled by Dufour (1984). One of the probable causes is our 
adopted temperature structure (equation ||) , which varies 
slightly from those used by others. For example, we find 
that our mean abundances for DEM L243 would increase 
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by about 0.1 dex if we adopted the relation of Campbell 
etal. (1986) at all values of T(0++). 

Our data are generally consistent with there being no 
abundance variations between the four different H n re- 
gions. It is interesting to note that DEM L199 is close to 
the LMC bar, about 1 kpc from the center of the galaxy; 
and DEM L243 is situated in the northern outskirts of the 
LMC-4 supergiant shell, at a galactocentric radius of ~3.5 
kpc. Pagel et al. (1978) have suggested that the LMC H n 
regions possibly exhibit a slight abundance gradient. This 
has not been further examined, nor has a gradient been de- 
tected in the cluster population (Olszewski et al. 1991). In 
our data, it is suggestive that DEM LI 99 and DEM L243 
delineate the extremes of any interpreted variation among 
our four objects. The difference in metallicity is consistent 
with the small gradient suggested by Pagel etal. (1978). 

3. SEMI-EMPIRICAL BRIGHT-LINE METHODS 

We turn now to examining more indirect emission-line 
diagnostics of metal abundances. In situations where T e 
cannot be adequately constrained by observation, it is 
common practice to estimate the metal abundances using 
the semi-empirical, "bright-line" abundance parameters. 
Here we examine the performance of these parameters in 
light of our detailed nebular data and highly-constrained, 
tailored photoionization models from Paper I. 

We also compute model tracks of the abundance param- 
eters, using Mappings II with generalized nebular param- 
eters. These incorporate the stellar energy distribution of 
CoStar model C2 (Schaerer & de Koter 1997), which cor- 
responds to an 06 - 07 stellar effective temperature. We 
assume an inner radius to the gas distribution of 0.4i?Sj 
and gas density n — 10cm -3 . In the figures that fol- 
low, the dashed, solid, and dotted lines correspond to the 
volume-averaged \ogU = —2, —3, and —4, respectively, 
which is equivalent to changing the total ionizing photon 
emission rate or gas filling factor. 

The grid of models is computed with Z = 0.05, 0.1, 
0.3, 0.5, 1.0, and 2.0 times Z Q . We included the el- 
ements (He, C, N, O, Ne, Mg, Al, Si, S, Ar, Ca, and 
Fe) with Z = (-1.01,-3.44,-3.95,-3.07,-3.91,-4.42, 
-5.53,-4.45,-4.79,-5.44,-5.88, and -4.96), respec- 
tively. We largely follow McGaugh (1991) in scaling the 
abundances of individual elements with respect to O. For 
He and N, we use the relations given by McGaugh, but 
scaled to match the Anders & Grevesse (1989) values for 
Z : 

He/H = 0.0850+ 15(0/H) (11) 

and 

log(N/H) = 1.5 log(0/H) + 0.66 . (12) 

For C and Fe, we adopt McGaugh's relations directly (his 
equations 10 and 11). The remainder of the elements are 
fixed in their proportion to O at Z , as given by Anders 
& Grevesse. Models for Z and 2 Z use the CoStar C2 
atmosphere at solar metallicity, while the rest use the cor- 
responding SMC metallicity model; we find that the stellar 
metallicity is unimportant for these ?7-tracks. 

3.1. i?23 

The O abundance parameter i?23 (Pagel etal. 1979; 
equation [j] above) , has been extensively used and empiri- 



cally calibrated several times, by McGaugh (1991), Skill- 
man (1989), and Dopita & Evans (1986), among others. In 
Figure 0a, we show log(0/H) vs. logi?23, with the tracks 
showing results from generalized photoionization models 
described above. Our assumptions differ somewhat from 
those of previous studies, in particular with the assump- 
tion of a fairly hollow morphology and new stellar atmo- 
sphere models. Most authors (e.g., McGaugh 1991; Dopita 
& Evans 1986) also assume some anticorrelation between 
the characteristic T+ or U, and Z, in their adopted cal- 
ibration at high abundance. Our tracks do not assume 
this anticorrelation. Despite these differences, our tracks 
remain similar to those of previous authors, although our 
models show a slight offset to higher log(0/H) (see Kob- 
ulnicky etal. 1999 and McGaugh 1991 for comparisons of 
i?23 calibrations). 

We also plot in Figure 0a the spatially integrated data 
for our objects, using the abundances derived from the 
means of our resolved apertures (Table ||) . The values of 
i?23 are computed in Paper I and shown here in Table 0. 
We caution that the nebular fractional area included in the 
spatial scans varies among the four objects, and we refer 
the reader to Paper I for the precise details. The three 
spatial scans of the spherical object DEM L323 (trian- 
gles) should give an indication of the degree to which sub- 
sampling is representative of the total spatial scan (solid 
triangle). The three measurements of i?23 are in excel- 
lent agreement, which is consistent with our finding in Pa- 
per I that this index is robust to spatial variations. For 
DEM L243, we show i?23 derived from the spectrum with 
the SNR-contaminated region subtracted (solid diamond) ; 
and also that from the total integrated region including 
the SNR (open diamond). The square and cross show 
DEM L199 and DEM L301, respectively. 

Our data points are generally well-behaved with respect 
to the model tracks in Figure 0a. While we found excellent 
agreement between the observed emission-line spectra and 
our tailored photoionization models in Paper I, wc sec in 
Figure 0a that most of the objects fall in their expected lo- 
cation with respect to the more generalized model tracks. 
DEM L323, DEM L243, and DEM L199 fall between tracks 
of log [7 = -3 and -4, with DEM L199 showing the high- 
est value of U, as expected in this high-excitation object. 
The one anomalous point is DEM L301 (cross), which we 
found in Paper I to have an unusual combination of exci- 
tation mechanisms. We concluded that this object, which 
has an extreme shell morphology, is most likely ionized 
by a combination of density-bounded photoionization plus 
shocks. While our tailored model for this object repro- 
duced the observed i?23 well, it is apparent that its value 
is anomalously high with respect to the tracks in Figure 0a. 
Ironically, the offset is in the sense of a higher ionization 
parameter, although the object in fact has a much lower U 
than the others in the sample. The larger value of i?23 is 
probably caused by the enhanced emission contributed by 
the shock activity. DEM L301 and DEM L199 have similar 
values of i?23, which, lacking any additional information 
for these objects, would imply similar abundances; we see 
that in fact this would overestimate log(0/H) for DEM 
L301 by about 0.3 dex, taking the measured log(0/H) at 
face value. 

We also overplot with small plus signs in Figure 0a the 
data compiled by Diaz & Perez-Montero (2000). These 
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Fig. 7. — Abundance parameters logi?23 (panel a) and log5"23 (panel b) vs. log(0/H) and log(S/H), respectively. The tracks indicate 
photoionization models for 0.05, 0.1, 0.3, 0.5, 1.0, and 2.0 Z©, as labeled in panel b. The dashed, solid, and dotted linetypes correspond to 
log U = —2, —3, and —4, respectively. Our observations are indicated by the symbols shown in the key; the solid triangle corresponds to 
observations for DEM L323 scanned over the entire object, while the open triangles represent partial scans (see Paper I). The small plus signs 
show data compiled by DPM, for objects whose abundances are determined using direct measurements of T e . 
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represent data from the literature (their Table 2) for which 
abundances were derived from a direct measurement of 
T[0 III]. Although the points for our data from Paper I 
are consistent with the model tracks, it is apparent that 
the photoionization models in Figure pp do not track well 
the locus of the larger dataset. It is important to note that 
adopting softer stellar atmospheres can improve the cor- 
respondence slightly, since this would offset the tracks to 
slightly lower log(0/H). We refer the reader to McGaugh 
(1991) to evaluate the consequences of the stellar effec- 
tive temperature. The atmospheres adopted here (CoStar 
C2) correspond to 06 - 07 stars, which are already cooler 
than most of our LMC objects, and only relatively small 
changes result if 03 - 04 atmospheres (CoStar E2) are 
used instead. The discrepancy between models and data 
has always been a difficulty in the use of R23, and therefore 
empirical calibrations of this parameter have been crucial 
for its successful use. 

3.2. ,523 

As mentioned earlier, a parameter similar to i?23 has 
been introduced for S by Vflchez & Esteban (1996) and 
Christensen etal. (1997), which was further explored by 
Diaz & Perez-Montero (2000, hereafter DPM). It is im- 
portant to note that this parameter, 523 (equation 0), is 
not strictly analogous to R23. While O and S have homol- 
ogous energy levels, the ionization potentials (IP) for their 
respective ions are different. In particular, it is important 
to note that the IP required to reach S +3 (34.83 eV) is 
virtually identical to that necessary for ++ (35.12 eV). 
Therefore, although S + and S ++ are indeed the dominant 
ions for S, for typical H n regions, there is likely to be non- 
negligible S + , which is ionized by the same radiation that 
produces ++ . Although Christensen et al. (1997) pointed 
out that the ionization fraction of S +3 is relatively small, 
typically <J 0.2, we show below that it nevertheless signif- 
icantly affects the ionization balance of S + and S ++ , and 
consequently, the value of 523. 

Figure 06 is similar to panel a, now showing log (S/H) 
vs log 523 (Table 0). The model line types and data sym- 
bols are the same as before. It is immediately apparent 
that, contrary to earlier claims in the literature, 523 is 
more sensitive to the ionization parameter than i?23. The 
change in log 523 between the model tracks, varying log U 
from -2 to -4, is almost 0.5 dex, whereas it is less than 0.3 
dex for i?23. The greater [/-sensitivity of 523 is caused by 
the "missing" contribution of S IV. Figure 06 shows that 
the models with high U show lower log 523, the opposite 
pattern to i?23. This is consistent with the ionization be- 
havior of S, since a larger population of S +3 is expected at 
higher U . 

As suggested by DPM, the lower- metallicity branch of 
523 does span a larger range in values than i?23, for a 
given range of Z (Figure 0). However, the location of the 
inflection at the maximum 523 is at only a slightly higher 
Z than that for i?23. Figure shows that our models for 
logi?23 have a maximum close to Z = 0.3 Z©, and those 
for log 523 have a maximum around Z = 0.5 Z©. Thus, 
there is only ~ 0.2 dex extension in the dynamic range of 
Z in the use of 523. Nevertheless, since so many of the 
observed objects in the literature have abundances around 
0.3 — 0.5 Z©, this augmentation makes a substantial differ- 
ence in evaluating abundances. As is dramatically shown 



in Figure 1 of DPM, 523 empirically shows an evidently 
monotonic increase as a function of Z, in contrast to R23, 
which shows a distinct double- valued structure in Z. 

Our data points are again reasonably consistent with 
the models in Figure 0, although they now perhaps show 
a tendency to fall between the \ogU = — 2 and —3 mod- 
els, rather than -3 to -4 for Figure 0a. The tailored 
photoionization models and spatially-resolved data in Pa- 
per I showed similar minor discrepancies. The data points 
for the shock-affected regions, DEM L301 (cross) and the 
SNR-contaminated observation for DEM L243 (open dia- 
mond) are offset to higher log 523. This suggests that 523 
is increased by the presence of shock excitation, similar to 
the behavior of R23 found in the previous section. 

We also see in Figure that the photoionization mod- 
els do track the data well for 523, and much better than 
for R23 in panel a. However, it is also apparent that the 
points with the highest values of 523 fall outside the model 
tracks. In formulating a calibration for 523, we would 
therefore recommend that these values be excluded, and 
that the empirical calibration derived by DPM should not 
be used for Z > 0.5 Z©. 

3.3. 5234 

As discussed in the previous section, the population of 
S +3 , which is ionized by the same radiation as that ion- 
izing ++ , is not sampled by the 523 parameter. While 
the ionization fraction of S +3 is relatively small, we have 
seen in the previous section that it significantly compro- 
mises the utility of 523 as an abundance diagnostic. We 
therefore suggest that the parameter, 

[S II] A6724 + [S III] AA9069, 9532 + [S IV] A10.5/Z 

(13) 

is a better abundance indicator than 523. In the same 
way that i?23 samples all significant ions of O, 5234 more 
completely samples the significant ions of S. Note that 
any population of S +4 (IP 47.30 eV) will be an insignifi- 
cant fraction of the total S ions, for massive star sources: 
S +4 /S <; 0.02 even for an object ionized by a WR star 
where He++/He = 0.4. The principal difficulty with 5234 
is the inclusion of the mid-IR line [S IV]A10.5/x, which is 
not readily observable with standard ground-based instru- 
mentation. However, since the IP necessary to produce 
S +3 is virtually the same as that for ++ , it is possible 
to estimate the abundance of S +3 based on that for ++ . 
This was demonstrated earlier by Mathis (1982) and Den- 
nefeld & Stasihska (1983). We present this approach here. 

Figure H presents log (S/H) vs. 5234, on the same scale 
as that for 523 in Figure 0b. The line types and symbols 
are the same as before. We see that 5234 is dramatically 
less sensitive to the ionization parameter, owing to the in- 
clusion of the S IV indicator. Indeed, the models for 5234 
are even less sensitive to U than is i?23 (Figure 0a), for 

Z< z©. 

It would therefore be desireable to estimate the in- 
tensity of [S IV] A 10. 5/i from that of the O ion- 
ization indicator, [O III]/[0 II]. Figure || shows 
log([S IV]A10.5/V[S III]AA9069,9532) vs. 
log([0 III]AA4959,5007/[O II]A3727) for models with an 
E2 CoStar atmosphere, and line types as before. This at- 
mosphere corresponds to an 03 - 04 stellar type, and we 
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Fig. 8. — Photoionization models and data for log(S/H) vs. S234. As before, the tracks are computed for metallicities 0.05, 0.1, 0.3, 0.5, 
1.0, and 2.0 Zq; and dashed, solid, and dotted linetypes correspond to logU = —2, —3, and —4, respectively. Observed values for 5234 are 
computed as described in the text, and are indicated by the symbols as in Figure M. 



prefer this model in examining the relation between S IV 
and other ions since it is more relevant in environments 
with harder ionizing fields. We see that for Z < 0.5 Z Q 
(solid points), the relation between these ratios is essen- 
tially a simple power law. For these points, we fit: 



log 



[SIV]A10.5/x 
[SIII]AA9069,9532 



-0.984+ 1.276 log 



[O II]A3727 
n (14) 

which is shown by the dot-dashed line in Figure O. 

With observations of [S III], this relation allows an es- 
timate of the [S IV] intensity, which can then be used 
to compute 5234. We note that the adoption of the 
cooler C2 CoStar atmospheres, used in our other pho- 
toionization models, would result in a difference of less 
than 0.02 and 0.01 in the fitted intercept and slope, re- 
spectively. Furthermore, the correction for [S IV] will 
only be significant for moderate to high U and/or high 
Z. As a test of equation we use mid-infrared and op- 
tical line observations of the Orion nebula by Lester, Din- 
erstein, fc Rank (1979). For their measurements, equa- 
tion |l4 predicts a volume emissivity for [S IV]A10.5/i 

-6 a — i which agrees within 



of 6.5 ± 2.7 x 10~ 21 erg cm - " S" 

measurement uncertainties with the observed value of 
9.0 ±2.7 x 10 -21 erg cm~ 6 s _1 . Considering the extremely 
narrow, 10" line of sight on the Orion nebula, and much 
higher density (10 4 cm~ 3 ) than considered for our pur- 
poses, this agreement is highly encouraging. 

It is thus relatively simple to convert from 523 into 5234 
and thereby almost eliminate the sensitivity to U. We used 
this method of estimating [S IV] to compute 5234 for our 
spatially- integrated observations (Table 0), which are plot- 



ted in Figure g, using the same symbols as before. In the 
errors for 5234, we include in quadrature an uncertainty 
of 25% for the uncertainty of [S IV] from equation |lj. We 
again have excellent agreement with the models. It is clear 
that Z can be estimated with greater confidence based on 
5234 than 523 at these metallicities, because the large 

OIII1AA4959 500fP reac ^ * n ^ e m °dels f° r ^23 (Figure |7p) has been vastly 
V ' ~" ' ' ' ' reduced for 5234. 

Similarly, we show in Figure [l^ that spatial variations 
are also reduced from 523 to 5234. Figures [Io| a and b show 
our spatially resolved observations of DEM L199 for 523 
and 5234. The solid line indicates the tailored model for 
this object, using the early WR model of Schmutz et al. 
(1992; see Paper I), central hole radius of 0.5i?s, an d gas 
density n = 100 cm~ 3 . While Figure ^0| a shows a large 
spatial variation of J> 0.4 dex for 523, we see in Figure [lOp 
that the variation in 5234 is reduced by about a factor 
of 2 in the logarithm. Figures |To| e — / show the same 
behavior for DEM L323. The solid line again represents 
the corresponding tailored model, with an 03 - 04 stellar 
atmosphere (Costar E2), central hole radius of 0.4i?g, and 
gas density n — 10 cm~ 3 . In Figure ^0| g — d, we show the 
spatially resolved data for DEM L199 superimposed on the 
model tracks for 523 and 5234, respectively. The reduced 
scatter in 5234 against the models again demonstrates the 
improved constraints in estimating log(S/H), compared to 
523. 

3.4. Calibrations 

In Figure [Tl|a and b we show the models for log(S/H) as 
a function of 523 and 5234, overplotted with the Galac- 
tic and LMC data presented by Dennefeld & Stasihska 



16 



Calibration of Abundance Diagnostics 



00 



> 

00 

en 

o 




1 .5 



1.0 -0.5 0.0 
log [0 lll]/[0 II] 



0.5 



1 .0 



Fig. 9.— Models of log([S IV]A10.5/V[S III]AA9069,9532) vs. log([0 III]AA4959,5007/[O II]A3727). These models employ the Costar E2 
stellar atmosphere model. The line types are as before, and the metallicity extremes are indicated in solar units for the U = —2 track (models 
for 0.05 and 0.1 Zq are virtually degenerate and indistinguishable). Equation (dot-dashed line) is fitted from the solid points (Z < 0.5 Zq). 



(1983). Their S abundances are computed from measure- 
ments of T[0 III] and observations of [S III]AA9069,9532. 
We compute 5234 from this dataset with the aid of equa- 
tion [l4|, as described above. Figure [TT] shows that the data 
present a well-defined sequence in both 523 and 5*234. It 
is especially encouraging that the locus of the models is 
in excellent agreement with that of the data, in contrast 
with the situation for i?23, as we saw above in Figure Ma. 
We replot the i?23 models with the Dennefeld & Stasinska 
data in Figure |Tl]c, again suggesting the same discrepancy 
seen earlier. It is apparent that for these data, the values 
of i?23 are fairly insensitive to log(0/H) as we saw before, 
owing to the location of the inflection and spread in U. 

To estimate a theoretical calibration for 523, we take 
the mean of the three models at each metallicity, up to 
0.5 Zq. A resulting power-law fit is shown by the lighter, 
straight, solid line in Figure |Tl]a. For 523, we obtain: 

log(S/H) = -5.43 + 1.33 log 523 . (15) 

The light dash-dot and dashed lines in Figure [Tl] a show the 
DPM and Christensen etal. (1997) calibrations, respec- 
tively. DPM calibrated a relation for log(0/H) vs. 523, 
so we used the solar S/O ratio to convert their relation 
to a calibration of log(S/H). It is apparent that all three 
calibrations are similar. The DPM relation shows the best 
correspondence to the data, as is expected since it is fitted 
to the largest dataset. It is especially encouraging that 
our theoretical relation is intermediate between the two 
empirical ones, confirming that the theoretical calibration 
is fully consistent with the available data. However, in us- 
ing any of these 523 calibrations, it is important to bear 



in mind that the models predict a double-valued relation 
around log 523 £ 0.0. 

We use the same procedure to fit a theoretical calibra- 
tion for 5234 and obtain: 

log(S/H) = -5.58+ 1.27 logS234 . (16) 

As in the case for 523, the data are in excellent agreement 
with this rough theoretical fit, shown by the light, solid 
line in Figure [Hp. 

We again emphasize that, although the data at Z > 
0.5 Zq are consistent with the calibrations for both 523 
and 5234, they strongly diverge from the models in that 
regime, and that a power-law approximation is necessarily 
crude near these values. We therefore consider the cali- 
brations reliable only for Z <J 0.5 Zq, and extreme caution 
should be exercised in extrapolating at higher metallicity. 
It is also essential to note that the double- valued nature of 
all of these abundance parameters still remains an issue. 

In Figure [Tl|d, we show the measured log(0/H) vs. 
log 5234 for the Dennefeld & Stasinska sample, where 
log(0/H) are again derived from direct measurements of 
T[0 III]. We see that the scatter is much larger than for 
log(S/H) vs. log 5234 (Figure |n]&). Although Garnett 
(1989), among others, suggests that there is no systematic 
variation in S/O with O/H, Figure [ll] shows that there is 
still significant variation in the S/O ratio among the dif- 
ferent objects. Therefore, while 5234 appears reasonably 
reliable for estimating the S abundance, it appears to be 
significantly less reliable for inferring the O abundance, 
and caution should be exercised accordingly. 
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Fig. 10.— Spatial variation of 523 and 5234 in DEM L199 (panels a and 6) and DEM L323 (panels e and /). The light, horizontal bars 
indicate the spatial extent of the apertures for these data. Panels c and d show the same data for DEM L199 superimposed on the models, 
with linetypes as in Figure [ 
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Fig. 11. — S abundance vs. diagnosticSjS'23 (panel a) and S23A (panel b), and O abundance vs. R23 (panel c) and S234 (panel d). The 
models are the same as in Figures |ij and 8l Data points are from Dennefeld & Stasinska (1983), with S234 computed as described in the 
text. The lighter, solid lines show a nt to our models (see text), while the light, dot-dashed and dashed lines show the DPM and Christensen 
etal. (1997) empirical calibrations, respectively. 
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4. CONCLUSION 

We have carried out a detailed investigation of elemen- 
tal abundance derivations using four H n regions in the 
LMC. We use tailored photoioinzation models to examine 
standard abundance analyses based on measured values of 
T e . Our data (Paper I) are derived from both spatially- 
resolved observations extracted from stationary long slit 
positions, and scanned, spatially-integrated slit observa- 
tions. We also examine the bright-line abundance diag- 
nostics for O and S, in light of the direct abundance de- 
terminations and photoionization models. 

Our abundance determinations are based on measure- 
ments of T[0 III], which we take to represent T(0 ++ ), 
and we assume a two-zone temperature structure for the 
nebulae, represented by T(0 ++ ) and T(0 + ). We use stan- 
dard ionic abundance relations to then determine the total 
elemental abundances for He, N, O, Ne, S, and Ar, with 
respect to H. Comparison with tailored Mappings pho- 
toionization models highlights the importance of choosing 
a relation between T(0 ++ ) and T(0 + ) that adequately 
represents the nebular temperature structure. Failure to 
do so can result in metallicity estimates that are discrepant 
by at least 0.2 dex from values indicated by tailored pho- 
toionization models. 

Abundance measurements for the stationary slit posi- 
tions show high spatial uniformity, with no evidence of 
variations or gradients to within 0.1 - 0.15 dex. Thus it 
is unlikely that there are systematic biases resulting from 
the strong ionization gradients seen in these objects. The 
adopted two-zone T e structure therefore appears to be 
highly reliable for estimating ionic abundance estimates 
even through "pencil-beam" apertures that sample only a 
small nebular area, at least for our fairly isothermal H n 
regions. 

No areas of local enrichment were detected in DEM 
L199, in spite of the presence of two WN3 stars and one 
WC4 star. The stellar products may be hidden in hot, 
coronal gas within the central superbubble, or the stars 
may not have produced enough enriched material to be 
readily detectable. The results show that self-enrichment 
by WR stars is likely to be a complex phenomenon, empir- 
ically. DEM L243 and DEM L301, both showing evidence 
of recent SNR activity, also do not show local enrichments, 
although with poorer constraints. 

Abundance measurements from the scanned, spatially 
integrated apertures are consistent with those obtained 
from the spatially resolved observations. Our results are 
~ 0.2 dex lower than average LMC Hll region measure- 
ments (Dufour 1984; Garnett 1999), probably resulting 
in part from different descriptions for the T e structure. 
The spatially-integrated measurements are also consis- 
tent with there being no variation between the four H n 
regions, although, interestingly, they are also consistent 
with the marginal abundance gradient suggested by Pagel 
etal. (1978). While the presence of the SNR in DEM L243 
did not affect the resulting abundances from the spatially 
integrated observation, the derived abundances for DEM 
L301 are on the low end of the distribution, hinting at 
spurious effects caused by the shock activity in that su- 
perbubble. 

We computed the i?23 O abundance parameters (Pagel 
et al. 1979) for the spatially integrated data, and compared 



these with model tracks constructed with Mappings. The 
models assume the Costar C2 (Schaerer & de Koter 1997) 
stellar atmosphere corresponding to an 06 - 07 spectral 
type, and an inner nebular radius of 0.4 Rg. As has histor- 
ically been the case, the models do not agree well with the 
locus of observations in the literature, although our LMC 
data do agree well, coincidentally, with both. 

Similarly, we examined the 523 abundance parameter 
for S (e.g., Christensen etal. 1997; Diaz & Perez-Montero 
2000). Our models reveal that, contrary to previous sug- 
gestions, 5*23 is more sensitive to the ionization parameter 
than is R23. S IV is produced by the same radiation that 
ionizes O III, and is a significant ion of S in many H n re- 
gions, but it is not sampled by 523. Its omission therefore 
causes 523 to be much more sensitive to U than R23. The 
spatially resolved observations confirm this by showing, in 
agreement with model predictions, lower values of 523 in 
the central nebular regions where S IV is important. As 
shown in Paper I, this spatial variation is not predicted or 
observed in i?23. 

Our models also suggest that the maximum in 523 oc- 
curs at only ~ 0.2 dex higher in Z than in R23. Neverthe- 
less, this appears to significantly alleviate the effect of the 
double- valued structure of log(S/H) vs 523 when inferring 
abundances, as shown by Diaz & Perez-Montero (2000). It 
is highly encouraging that the data, both from our sample 
and from the literature, are in excellent agreement with 
the models, in contrast to the behavior of i?23. We offer a 
theoretical calibration of 523 (equation |l5|) which appears 
to be fully compatible with the data in the literature thus 
far. However, we caution that the locus of the available 
data may well be deceptive in suggesting that a power-law 
relation can be used at Z ;> 0.5 Z Q . 

To overcome the limitations of 523 in [/-sensitivity and 
spatial variation, we introduce a similar S abundance pa- 
rameter, 5234. This is the same as 523 with the added 
emission of [S IV]A10.5//. Although this mid-IR line is not 
readily observable with most conventional ground-based 
spectrographs, it is straightforward to estimate its inten- 
sity from the simple correspondence between [S IV] /[S III] 
and [O HI]/[0 II] (equation |lj). Our models show that 
5234 is less dependent on U than is even i?23. 5234 for our 
objects and for the larger sample of Dennefeld & Stasihska 
(1983) are in excellent agreement with the models. Like- 
wise, the spatial variations for both models and observa- 
tions are dramatically reduced for 5234 in contrast to 523. 
We provide a theoretical calibration for log(S/H) vs 5234 
at Z < 0.5 Zq (equation 16). 

Finally, we reiterate some caveats for the use of i?23, 
523, and 5234. We find that the presence of shock ex- 
citation increases the value of these parameters; for our 
objects, the effect is about 0.1 dex in magnitude. Sec- 
ondly, significant variations in the S/O ratio dictate cau- 
tion in inferring O abundances using 5234 and 523 (Fig- 
ure |TT| ) . It is also important to bear in mind the double- 
valued structure for all three of these parameters. Lastly, 
we emphasize the deviation between the data and models 
above 0.5 Z , and we therefore consider the calibrations 
presented thus far for 523 and 5234 to be reliable only 
for Z <J 0.5 Z . Further empirical investigation is needed 
to understand the behavior of these parameters at higher 
metallicity. Bearing in mind these caveats, the excellent 
correspondence between the modeled 5234, 523, and the 
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available data, together with the more highly monotonic 
behavior of these parameters, promises greater effective- 
ness as metallicity indicators than R23. With improving 
access to the [S IV]A10.5^ line, it should be possible to 
confirm the behavior of 5234 directly. 

It is a pleasure to acknowledge discussions with Mike 
Dopita, Annette Ferguson, Don Garnett, Dick Shaw, Evan 
Skillman, Elena Terlevich, and Bob Williams. We are also 
grateful to Mike Dopita for access to the Mappings II 
photoionization code and to Angelez Diaz for access to her 
work in advance of publication. Finally, we are pleased to 
acknowledge the referee, Bernard Pagel. 
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Table 1 

Hll REGION SAMPLE 



DEM 


Henize 


LH 


Sp. Type 


log(L Ha /crg s x ) a 


Inner radius 


Shocks 


Sp. Type Reference 


DEM L199 


N144 


LH 58 


WN3 


38.6 


0.5 i? s b 


no 


Breysacher (1981) 


DEM L243 


N63 A 


LH 83 


07 


37.7 


0.1 R s 


yes 


Oey (1996a) 


DEM L301 


N70 


LH 114 


03 


37. 7 C 


0.95 R s 


yes 


Oey (1996a) 


DEM L323 


N180 B 


LH 117 


03-4 


38. 4 d 


0.4 R s 


no 


Massey etal. (1989) 



a From Oey & Kennicutt (1997), unless otherwise indicated. 
b DEM L199 has complex morphology (see Paper I). 
C DEM L301 is significantly density-bounded (Paper I). 
d Kcnnicutt (1997), private communication. 
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Table 2 

Logarithmic abundances of spatially resolved observations' 1 



Observation [Ho/H] b orr [N/H] err [O/H] err [Ne/H] c err [S/H] err [Ar/H] err 



D199.205-apl0 


-1 


.059 





026 


4 


.892 





.126 


3 


.647 





.216 


-4. 


390 





256 


5 


.282 


0. 


116 


-5 


.993 


0. 


133 


D199.205-apll 


-1 


.091 





.021 


-4 


899 


0. 


.044 


-3 


.633 





.065 


-4. 


374 





.076 


-5 


.261 


0. 


038 


-5 


.896 


0. 


046 


D199.205-apl3 


-1 


.039 


0. 


.015 


-4 


.890 





.040 


-3 


.621 





.056 


-4. 


368 





.064 


-5 


.280 


0. 


032 


-5 


.889 


0. 


040 


D199.205-apl4 


-1 


.066 


0. 


020 


-4 


.977 


0. 


.076 


-3 


571 





130 


-4. 


.409 





158 


-5 


.172 


{) 


063 


-5 


.583 


0. 


080 


D199.205-apl5 


1 


.057 





026 


-4 


.927 


0. 


.118 


-3 


.681 


0. 


204 


-4. 


483 





241 


-5 


.378 


0. 


105 


-5 


.886 


0. 


125 


D199.205-apl6 


-1 


.016 


0. 


.027 


-4 


.946 





175 


-3 


.421 


0. 


.313 


4 


080 


{) 


372 


-5 


.315 


0. 


157 


-5 


.790 


0. 


186 


D199.205-apl7 


-1 


.081 


0. 


.029 


-4 


.947 


0. 


.104 


-3 


.689 


0. 


.182 


-4. 


.537 


0. 


225 


-5 


.315 


0. 


096 


-5 


.924 


0. 


113 


D199.205-apl8 


-1 


.001 





069 


-4 


.545 





.327 


-3 


.018 





.587 


-3. 


.635 


0. 


706 


-4 


.958 


0. 


.308 


-5 


.577 


0. 


358 


D199.205N120-aplO 


-1 


.065 


0. 


029 


-4. 


.928 


0. 


.178 


-3 


.669 


0. 


.302 


-4 


.316 


0. 


.359 


-5 


.316 


0. 


.161 


-6. 


.010 


0. 


186 


D199.205N120-apll 


-1 


.065 





027 


-4. 


.932 


0. 


.137 


-3 


.627 


{) 


.237 


-4 


.337 


0. 


.283 


-5 


.337 


0. 


125 


-5. 


.960 


0. 


147 


D199.205N120-apl2 


1 


.083 


0. 


015 


-4 


.999 


0. 


043 


-3 


771 


0. 


063 


-4 


.553 





074 


-5 


.391 


0. 


037 


-6. 


.104 


0. 


047 


D199.205N120-apl3 


-1 


.088 


0. 


018 


-4 


.950 


0. 


046 


-3 


695 





072 


-4 


.423 


0. 


085 


-5 


.337 





041 


-5 


.992 


0. 


049 


D199.205N120-apl4 


-1 


.011 


0. 


015 


-4 


.898 


0. 


.047 


-3 


.552 


0. 


.073 


-4 


.317 





.085 


-5 


.268 


0. 


040 


-5. 


.814 


0. 


049 


D199.205N120-apl5 


-1 


.033 


0. 


019 


-4. 


.942 





.077 


-3 


.573 


{) 


.130 


-4 


.344 


{) 


.154 


-5 


.316 


0. 


.068 


-5. 


.867 


0. 


081 


D199.205N120-apl6 


-1 


.031 


0. 


019 


-4 


.869 





.089 


-3 


.569 





154 


-4 


.347 





184 


-5 


.259 


{) 


.081 


-5 


.815 


0. 


097 


D199.205N120-apl7 


-1 


.033 


0. 


027 


-5. 


.205 


0. 


.161 


-3 


.522 


0. 


296 


--4 


.271 


0. 


.334 


-5 


.264 


0. 


118 


-5 


.674 


0. 


148 


D199.205N120-apl8 


1 


.046 


0. 


027 


-4. 


.905 


0. 


141 


-3 


.382 





252 


-4 


.071 


{) 


.301 


-5 


.258 


0. 


.128 


-5. 


.755 


0. 


151 


D199.205N120-apl9 


1 


.014 


0. 


024 


-4. 


.909 





.131 


-3 


578 





230 


-4 


.304 





275 


-5 


.218 


0. 


119 


-5. 


.763 


0. 


143 


D243.2S-ap6 


1 


.068 


0. 


026 


-5 


.412 





.110 


-4 


177 





.175 


-5. 


.306 





242 


-5 


.710 


0. 


095 


-6 


.360 


0. 


115 


D243.2S-ap9 


-1 


.176 


{) 


.025 


-5 


.279 


0. 


151 


-4 


.094 


0. 


.257 


-5. 


.064 





.303 


-5 


.683 


0. 


136 


-6 


.393 


0. 


158 


D243.2S-apl4 


-1 


.313 


0. 


.031 


-5. 


.252 


0. 


.108 


-4. 


.055 


0. 


.173 


4. 


940 


0. 


.206 


-5 


.643 


0. 


.098 


-6 


.570 


0. 


119 


D243.5S-ap3 


1 


.559 


0. 


069 


-5. 


.408 


0. 


.350 


-4. 


.337 





563 


-5. 


442 





.775 


-5 


.801 


0. 


317 


-7 


.030 


0. 


338 


D243.5S-ap6 


-1 


.296 


0. 


023 


-5 


.079 





.177 


-3 


703 





305 


-4. 


675 





.363 


-5 


.414 





.159 


-6 


.200 


0. 


182 


D243.5S-ap5 


-1 


.101 





.019 


-5. 


.074 


0. 


.064 


-3 


.587 


0. 


.107 


-4. 


457 





.130 


-5 


.375 


0. 


.056 


-5 


.931 


0. 


066 


D243.5S-ap7 


-1 


.156 


{) 


.034 


-4. 


.843 





.290 


-3 


.315 





.507 


-4. 


215 





.607 


-5 


.246 


0. 


.262 


-5 


.885 


0. 


297 


D243.5S-ap8 


-1 


161 


0. 


.039 


-5 


.085 





.263 


-3 


634 





452 


4. 


.654 





.556 


-5 


.406 





235 


-6 


.073 


0. 


276 


D243.30S-ap9 


-1 


.293 





.125 


-5. 


.412 


0. 


.319 


-4. 


.397 


0. 


.520 


-5. 


.678 


0. 


.720 


-5 


.803 


0. 


298 


-6 


.617 


0. 


323 


D243.30S-aplO 


-1 


.146 


{) 


.047 


-5 


.219 





189 


-3 


.944 


0. 


.317 


-4. 


951 


{) 


.400 


-5 


.526 


0. 


169 


-6 


.242 


0. 


198 


D243.30S-apll 


-1 


.164 





036 


-4 


.763 





.298 


-3 


147 





.522 


-4. 


.079 





641 


-5 


191 


{) 


270 


-5 


.769 


0. 


312 


D243.30S-apl4 


1 


.257 





077 


-5. 


.413 





.240 


-4 


332 


0. 


.397 


-5. 


398 





524 


-5 


.786 


0. 


222 


-6 


.594 


0. 


254 


D243.30S-apl5 


-1 


.108 


{) 


.098 


-4. 


.587 


0. 


.320 


-2. 


.909 





.549 


-3. 


.563 


{) 


.675 


-4 


.980 


0. 


.292 


-5 


.608 


0. 


340 


D243fix.2S-ap6 


-1 


.083 





.025 


-5. 


.145 


0. 


.065 


-3 


.745 


0. 


.100 


-4. 


.796 





.150 


-5 


.476 


0. 


054 


-6 


.072 


0. 


068 


D243fix.2S-ap9 


1 


.185 


0. 


024 


-5. 


.083 





055 


-3 


.778 


0. 


089 


-4. 


687 


0. 


.116 


-5 


.508 


0. 


.048 


-6 


.180 


0. 


058 


D243fix.2S-apl4 


-1 


.332 





030 


-4 


.991 


0. 


056 


-3 


.633 





091 


-4. 


.441 





.114 


-5 


.405 


0. 


.051 


-6 


.289 


0. 


071 


D243fix.5S-ap3 


-1 


.575 





.079 


-5. 


.085 


0. 


.057 


-3 


.813 


0. 


.091 


-4. 


.824 


0. 


.411 


-5 


.513 


0. 


.048 


-6 


.693 


0. 


095 


D243fix.5S-ap6 


-1 


.295 


0. 


.022 


-5. 


.095 


0. 


.054 


-3 


.731 


0. 


.089 


-4. 


708 


{) 


.109 


-5 


.428 


0. 


049 


-6 


.216 


0. 


058 


D243fix.5S-ap5 


-1 


.096 


0. 


019 


-5 


.167 


0. 


.055 


-3 


.755 


0. 


.087 


4. 


.657 





109 


-5 


.459 


{) 


049 


-6 


.031 


0. 


058 


D243fix.5S-ap7 


-1 


.143 


{) 


.031 


-5. 


.096 


0. 


056 


-3 


.757 


0. 


090 


-4. 


.737 





.135 


-5 


.473 


0. 


.053 


-6 


.144 


0. 


066 


D243fix.5S-ap8 


1 


.158 


0. 


.037 


-5 


.130 


0. 


059 


-3 


.713 


0. 


089 


4 


747 


{) 


160 


-5 


.447 


0. 


052 


-6 


.120 


0. 


075 


D243fix.30S-ap9 


1 


.305 


0. 


132 


-5 


.010 





.058 


-3 


748 





.092 


-4. 


919 





432 


-5 


.446 


{) 


.055 


-6 


.196 


0. 


090 


D243fix.30S-aplO 


1 


.150 





047 


-5. 


.078 





.056 


-3 


715 





089 


-4. 


678 


0. 


.176 


-5 


.400 


0. 


050 


-6. 


.089 


0. 


063 


D243fix.30S-apll 


-1 


.151 


{) 


.033 


-5. 


.086 


0. 


.054 


-3 


.722 


0. 


.087 


-4. 


760 





.192 


-5 


.484 


0. 


050 


-6. 


.107 


0. 


060 


D243fix.30S-apl4 


1 


.283 


0. 


.076 


-5 


.089 





059 


-3 


.810 


0. 


.092 


-4. 


.783 


0. 


.278 


-5 


.493 


0. 


.059 


-6. 


.247 


0. 


075 


Uz4torix.oUo-a.pio 


1 


.084 





097 


—5. 


063 





065 


—3 


759 





093 


A 

—4. 


573 





236 


—5 


,42 1 


f). 


059 


—6 


105 


0. 


092 


D301.SW6-ap7 


-1 


.061 





019 


-5 


193 


0. 


059 


-3 


984 


0. 


.092 


-4. 


.695 


0. 


111 


-5 


.487 





052 


-6 


.323 


0. 


063 


D301.SW6-ap8 


-1 


.031 


0. 


.020 


-5 


.116 


0. 


.062 


-3 


.874 


0. 


.097 


-4. 


672 


0. 


.118 


-5 


.445 


0. 


055 


-6 


.205 


0. 


067 


D301.SW6-ap9 


-1 


.065 


0. 


020 


-4 


.975 


0. 


.199 


-3 


667 





336 


-4. 


377 


0. 


397 


-5 


.279 


0. 


178 


-6 


.052 


0. 


199 


D301.SW6-aplO 


-0 


.920 


0. 


.025 


-5 


.071 


0. 


.064 


-3 


.820 


0. 


101 


-4. 


475 





120 


-5 


.411 





057 


-6 


144 


0. 


070 


D301.SW6-apll 


-1 


.036 





.020 


-5. 


.163 


0. 


062 


-3 


.912 


0. 


.098 


-4. 


.743 


{) 


.125 


-5 


.439 


0. 


055 


-6 


.228 


0. 


067 


D301.SWl-ap6 


1 


.135 


{) 


.032 


-5 


.008 


0. 


309 


-3 


.710 


0. 


.514 


-4. 


190 


0. 


.607 


-5 


.362 


0. 


277 


-6 


.254 


0. 


305 


D301.SWl-ap5 


-1 


.088 





018 


-4 


.980 





195 


-3 


518 





334 


-4. 


437 


0. 


395 


-5 


.297 


{) 


.169 


-5 


.984 


0. 


193 


D323.C2-ap7 


1 


.194 





022 


-5 


.127 





.105 


-3 


.597 





.180 


-4. 


519 





213 


-5 


.327 


0. 


.090 


-5 


.958 


0. 


109 


D323.C2-ap8 


-1 


.089 


{) 


.017 


-5 


.168 


0. 


.070 


-3 


.710 





.114 


-4. 


625 





.135 


-5 


.421 


0. 


060 


-6 


.040 


0. 


074 


D323.C2-ap9 


-1 


.072 


0. 


.016 


-5. 


.084 


0. 


.048 


-3 


.594 


0. 


.078 


-4. 


437 





.091 


-5 


.349 


0. 


042 


-5 


.942 


0. 


051 


D323.C2-aplO 


1 


.190 





021 


-5. 


.113 


0. 


042 


-3 


.627 





060 


-4. 


483 


0. 


069 


-5 


.323 


0. 


.035 


-5 


.908 


0. 


042 


D323.C2-apll 


1 


.049 


{) 


.016 


-5 


.149 


0. 


.058 


-3 


.639 





.089 


-4. 


488 


0. 


.102 


-5 


.316 


0. 


.046 


-5 


.843 


0. 


056 


D323.C2-apl2 


-1 


.057 


0. 


.016 


-5 


.127 


0. 


.053 


-3 


.631 


0. 


082 


4. 


.511 


0. 


.096 


-5 


.294 


0. 


.043 


-5 


.883 


0. 


051 


D323.C2-apl3 


1 


.066 





027 


-5. 


.268 


0. 


083 


-3 


852 


0. 


.130 


-4. 


782 


0. 


159 


-5 


.419 


0. 


.072 


-6 


.074 


0. 


088 


D323.Cl-apl3 


-1 


.080 


0. 


.051 


4 


.995 





.302 


-3 


.461 





.526 


-4. 


.341 


0. 


639 


-5 


.265 


{) 


262 


-5 


.885 


0. 


313 


D323.Cl-ap6 


-1 


.169 


{) 


.019 


-4. 


.919 


0. 


.127 


-3 


.361 


0. 


.219 


-4. 


255 





.259 


-5 


.207 


0. 


108 


-5 


.912 


0. 


127 


D323.Cl-ap7 


-1 


.175 


{) 


.022 


-5 


.238 


0. 


053 


-3 


.660 


0. 


.088 


-4. 


.575 





.104 


-5 


.345 


0. 


.046 


-5 


.891 


0. 


054 


D323.Cl-ap8 


-1 


.050 





016 


-5 


.178 





052 


-3 


635 





082 


-4. 


513 





095 


-5 


.306 





.043 


-5 


.835 


0. 


051 


D323.Cl-ap9 


1 


.077 





016 


-5 


.165 





.033 


-3 


602 





.044 


-4. 


454 





.050 


-5 


.361 


0. 


028 


-5 


.920 


0. 


032 


D323.Cl-aplO 


-1 


.085 


{) 


.017 


-5 


.112 


0. 


.067 


-3 


.618 





.112 


-4. 


.525 


0. 


.132 


-5 


.359 


0. 


.058 


-5 


.993 


0. 


070 


D323.Cl-apll 


-1 


.059 





.016 


-5. 


.133 


0. 


.040 


-3 


.678 


0. 


.063 


4. 


.763 


0. 


.080 


-5 


.316 


0. 


036 


-5 


.995 


0. 


042 


D323.Cl-apl2 


-1 


.113 


0. 


029 


-5. 


.149 





253 


-3 


713 


0. 


425 


-4. 


530 


0. 


513 


-5 


.391 


0. 


217 


-6 


.159 


0. 


252 



a D243fix observations have values computed for T(0 + + ) fixed at 9700 K. 

b [He/H] are lower limits for DEM L243. 

c [Ne/H] is subject to systematic uncertainties (see text). 
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Table 3 

Mean logarithmic abundances from spatially resolved observations 



Object 


T[0 III](K) 


cr 


[Hc/H] 


<j 


[N/H] 




[O/H] 


<7 


[Ne/H] a 


<j 


[S/H] 


cr 


[Ar/H] 


a 


DEM L199 


9620 


140 


-1.05 


0.005 


-4.93 


0.017 


-3.64 


0.026 


-4.40 


0.030 


-5.30 


0.014 


-5.90 


0.018 


DEM L243 


9620 


360 


-1.18 b 


0.009 


-5.17 


0.040 


-3.83 


0.068 


-4.75 


0.084 


-5.50 


0.036 


-6.18 


0.043 


DEM L243 c 


9700 




-1.18 b 


0.009 


-5.09 


0.016 


-3.74 


0.025 


-4.67 


0.042 


-5.46 


0.014 


-6.17 


0.019 


DEM L301 


11900 


420 


-1.05 


0.008 


-5.13 


0.030 


-3.89 


0.047 


-4.63 


0.058 


-5.44 


0.027 


-6.22 


0.032 


DEM L323 


9480 


110 


-1.09 


0.005 


-5.15 


0.015 


-3.63 


0.022 


-4.52 


0.026 


-5.33 


0.012 


-5.93 


0.015 



a [Nc/H] is subject to systematic uncertainties (sec text). 
b [Hc/H] is a lower limit in this object. 
C T(0+ + ) fixed at 9700 K for all apertures. 



Table 4 

Logarithmic abundances from spatially integrated observations 



Observation 


T[0 III](K) 


err 


[Hc/H] 


err 


[N/H] 


err 


[O/H] 


err 


[Nc/H] a 


err 


[S/H] 


err 


[Ar/H] 


err 


D199.496W240 


8100 


1100 


-1.13 


0.03 


-4.77 


0.17 


-3.31 


0.31 


-3.97 


0.37 


-5.19 


0.16 


-5.73 


0.18 


D243.2(total) 


11900 


1000 


-1.18 


0.02 


-5.24 


0.08 


-4.01 


0.12 


-4.90 


0.14 


-5.54 


0.07 


-6.36 


0.08 


D243.2(no SNR) 


11400 


1400 


-1.17 


0.03 


-5.26 


0.12 


-3.97 


0.19 


-4.90 


0.23 


-5.58 


0.10 


-6.26 


0.12 


D301.SW6 


13000 


900 


-1.04 


0.02 


-5.19 


0.06 


-3.94 


0.09 


-4.66 


0.10 


-5.50 


0.05 


-6.27 


0.06 


D323.140(total) 


9600 


700 


-1.09 


0.02 


-5.10 


0.08 


-3.64 


0.14 


-4.54 


0.17 


-5.35 


0.07 


-5.97 


0.09 


D323.140N30N 


9700 


1400 


-1.11 


0.03 


-5.13 


0.18 


-3.66 


0.30 


-4.57 


0.36 


-5.35 


0.15 


-6.04 


0.18 


D323.140N30S 


8800 


400 


-1.21 


0.02 


-5.02 


0.06 


-3.48 


0.10 


-4.30 


0.12 


-5.27 


0.05 


-5.85 


0.06 


Dufour (1984) b 






-1.07 


0.02 


-5.03 


0.10 


-3.57 


0.08 


-4.36 


0.10 


-5.15 


0.11 


-5.80 


0.06 



a [Ne/H] is subject to systematic uncertainties (sec text). 

b Mcan LMC Hll region abundances as compiled by Dufour (1984). 



Table 5 

Abundance parameters for spatially integrated observations 



Observation 


i?,23 a 


err 


S23 a 


err 


S234 b 


err 


D199.496W240 


6.9 


0.4 


1.11 


0.08 


1.38 


0.38 


D243.2 (total) 


4.7 


0.5 


1.31 


0.08 


1.35 


0.36 


D243.2 (no SNR) 


4.7 


0.5 


1.04 


0.07 


1.09 


0.29 


D301.SW6 


7.1 


0.9 


1.83 


0.11 


1.85 


0.49 


D323.140 


6.1 


0.5 


1.19 


0.08 


1.31 


0.34 


D323.140N30N 


5.8 


0.6 


1.27 


0.09 


1.34 


0.35 


D323.140N30S 


6.2 


0.5 


1.15 


0.08 


1.29 


0.34 



a From Paper I. 

b Computed by estimating [S IV] with equation [l 
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